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Abstract 



> 

in 

i Cosmic rays of energy higher than lO^^el/ may be explained by 

00 topological defects produced in the early stages of universe. Two 

suitable alternatives are: necklaces formed by magnetic monopoles 



, connected by strings, and vortons which are loops stabilized by super- 



conducting currents. The former are uniformly distributed in the uni- 
Q^' verse, may account for cosmic rays above the ankle, suffer a transient 

Q . GZK cutoff with a subsequent recovery and isotropy of observations is 

I expected. The latter are concentrated in the galactic halo, require an 

^ ■ additional extragalactic contribution between the ankle and the GZK 

cutoff, beyond which give a harder component and predict anisotropy 



related to mass concentration. 
PACS: 98.70.S , 98.80. C 



1 Introduction 

The ultra-high energy cosmic rays (UHECR) have an energy spectrum of 
their flux F{E) that shows for E^F{E) a minimum around 5 x lO^^^el^ which 
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is called the ankle, then a maximum before the GZK cutoff|]T[ at 5 x lO^^eV 
due to the interaction with the CBR and a recovery after it at lO^^eV. 

Whereas the cosmic rays below the ankle are most probably of galactic 
origin, it is not clear which is the explanation of the subsequent rise. The 
possibihty that it is due to an extragalactic source seems to be supported by 
a partial GZK cutoff which would affect cosmic rays traveling at least 50Mpc. 
But the subsequent observed spectrum up to the highest energy event of 3 
X lO'^^eV indicates a hard component which may be or not related to that 
which appears above the ankle. 

It is difficult to explain the observed events beyond the GZK cutoff 
with ordinary astrophysical objects which are not identified if they are close 
to us, and that would require non standard messengers not interacting with 
CBR if they are very far away[0 . 

A solution may be the top-down mechanism where some type of super- 
heavy microscopical object with mass of the order of the Grand Unification 
Theory (GUT) scale decays very slowly producing the observed UHECR. 
Since a general feature is a hard spectrum at emission F(i?em) oc 1/Eem , 
these superheavy objects might be either condensed in the galactic halo or 
uniformly distributed in the universe. 

The former case may correspond to closed cosmic strings stabilized by 
superconducting currents called vortons|^ or to superheavy particles whose 
interaction with the ordinary ones is of gravitational order generically de- 
noted as cryptons|^ . Since they were presumably produced at the GUT 
scale, they behave afterwards as cold dark matter (CDM) and should have 
concentrated in the galactic halo. 

The latter case may be instead represented by necklaces where ordinary 
cosmic strings whose dynamics makes them evolve to a scaling solution, i.e. 
a uniform distribution in space, incorporate monopoles and antimonopoles 
that annihilate very slowly. 

We will give a brief description of vortons and particularly show that 
if they represent a small fraction of the halo CDM, they may account for 
the apparent hard component of the UHECR spectrum above lO^^eV. An 
additional extragalactic component would be necessary to explain the ankle 
feature. 

Compared to this, an also small fraction of the critical density of universe 
represented by necklaces may give, with a reasonable law for the energy de- 
grading due to interaction with CBR, the maximum of E^F{E) immediately 
below the GZK cutoff. The hard spectrum at emission allows a recovery 
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above it at least up to the highest observed UHECR which would be im- 
possible for extragalactic sources with the ordinary law F{Eem) oc l/E^^ 



2 Vortons in halo 

Considering sources that emit hx (t) GUT boson particles X per unit 
space and time, each of them giving Nc UHECR, the total flux on earth will 
be 
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F=±f\tN^n,it)(m . (1) 



where a is the scale parameter of universe, to its age and tin the initial time 
of contributions. 

For quasistable objects like vortons concentrated in halo with density n{t) 
which emit by tunneling an X with a lifetime r, nx= n/r and the total 
flux is 

^ Nc , ,At 

F^ = -^n^,{to)— , (2) 

where At ~ 50kpc is the halo size. 
The energy spectrum 

1 At 

F,{E) = n,{to)—Y.6{E-E,) (3) 
Att t .^^ 

must be averaged on the intervals AEi between produced particles to compare 
with observations 

F.{E.) = ^ MM ^ . (4) 
AE, An T ^ ' 

Since the production of UHECR comes from the hadronization of the very 
energetic quark into which the X particle decays, by dimensional arguments 
there are no relevant ordinary mass parameters and we may expect 

AE, ~ E, , Fh{E,) oc ^ , (5) 

i.e. a hard spectrum consistent with accurate QCD calculations 0] apart from 
energies close to mx- 
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It is interesting that according to Eq.(5) the average probabihty of UHECR 
production on energy intervals will be 



^^11 (6) 
dE T E ' 

which integrated on the ultra-high energy range 10^^ — lO'^'^eV to agree with 
the total flux Eq.(2) must give 



with A^'c ~ 10 , that is a reasonably accepted value [§] . 

Therefore we may take the equally spaced particles in log E according to 
El ~ lO^^eV, E2 ~ W^-^eV, E^ ^ W^eV, E^ ~ W^-^eV. Ei^ ~ 

1023.5g-j/ 

SO that AEi ~ 0.7£^i . 

Since from eq.(4) the flux in each energy bin is the same, one may nor- 
malize it roughly at the expected value for IQ'^^eV 

nhito) At ^ 1 
Arc r km? century 

A vorton is a loop of ordinary cosmic string with an energy per unit length 
/i ~ m\ stabilized by massless fermionic carriers giving therefore a total 
energy 

E, = fiL + —, (9) 

which is minimized at E^ ~ 2Nmx by a length L ~ N/mx ■ The decay 
of the vorton with emission of X by tunneling gives a lifetime r ~ to for 

~ 1000. Therefore to satisfy Eq.(8) one only needs a fraction ~ 10~^ of 
the average energy density of the halo 0.3 GeV/cm? represented by vortons. 
This small density may be the remnant of the collapse of most vortons at 
the electroweak phase transition ^ . 

This contribution of vortons to E^F{E) will be a hard component which 
reproduces the observed flux above the GZK cutoff as is seen in fig.l. One 
must therefore complement it with another possibly extragalactic component 
which may explain the spectrum between the ankle and GZK energy. 

A similar analysis may be done for cryptons being necessary to explain 
their required density. 
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3 Necklaces in universe 



These hybrid topological defects may be formed by a sequence of GUT sym- 
metry breakings 

G^Hx U(l) ^HxZ2 , (10) 

where in the first monopoles would be produced and then would be attached 
to the ordinary strings which appear in the second one as beads of a necklace. 

The relevant parameter for the necklace dynamics is r = m/{fid) where 
m is the monopole mass, d its separation from the antimonopole in the string 
and iJ, the tension of the latter. For r ~ 10® the distance between strings at 
present is small ~ SMpc. 

The evolution of the necklace networks is scale invariant, i.e. they would 
be distributed uniformly in the universe and represent a constant fraction 
of its energy. Monopoles and antimonopoles trapped in the necklaces at the 
end would annihilate producing X particles with a rate 

nx (t) ~ = -J • 11 

Therefore the expression for the UHECR flux eq.(l) would apply but 
with an early tj„ compatible with avoiding their redshift below lO^^eV which 
is roughly of the order of the matter-radiation equivalence time, r^/x can- 
not be larger than lO^^GeF^ to avoid a diffuse gamma radiation above the 
experimental bound. The total flux will be 

F„ = iV " In (^) . (12) 

Even though at emission the UHECR produced by an X are equally 
spaced in logE^g^ as discussed in Sec. 2, their redshift would cause a softening 
of the law XjE for the flux spectrum on earth. 

But more important than this effect is the interaction of cosmic rays with 
CBR. The p7 total cross-section at the highest energy is ~ 0.2 rah , and rises 
up to ~ 0.6 mh for the A resonance mass. 

Then to evaluate the flux spectrum, we may proceed as follows. Instead 
of taking one cosmic ray in each bin as in the case of vortons in halo, we will 
consider to account for the degrading of energy so that 
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Therefore we may parameterize for all cases 



log EfF{E,) = log J + log N, + 2 log 



( 



) 



(14) 



where J , related to the properties of sources, will be adjusted to fit the 
observed events. For the case of vortons all Ni = 1. 

For necklaces to determine the effective Ni and considering the mean free 
path associated to the quoted cxp^ and the density of 3K photons, we will 
take 1% of probability that the cosmic ray for A production keeps its energy. 
This is consistent with the fact that the sources of non degrading protons are 
concentrated in a radius ~ 50Mpc out a whole space ~ 100 times larger. For 
higher Ei this probability will increase to 3% following ap-y. We will assume 
that the missing events are transferred in equal parts to the two inmediate 
lower bins. According to the bins defined in Sec. 2, that of E^ is particularly 
affected by the resonant scattering. The inmediate lower one has only the 
upper 10% of the bin in the resonance region so that ~ 90% of its events 
keep their energy and the rest is transferred to the bin of Ei . 

In this way it turns out 

Ni = 4.09 N2 = 5.15 N3 = 0.05 = 0.14 iVg = 0.12 
Ne = 0.10 Nj = 0.08 iVg = 0.06 A^'g = 0.04 = 0.03 
and the flux spectrum eq.(14) is shown in fig.2. 

We see that the existence of an ankle and the recovery after a transient 
GZK are successfully reproduced. 

A check of the calculation is that the normalization for the bin of lO^^el^ 
must be 



which is satisfied for a = 10^''sec~^ coming from mx = lO^^GeV. 

This is similar to the normalization for ordinary strings |P the difference 
being that for them /x ~ mj^ and the present separation between strings is 
three orders of magnitude larger than for necklaces. As a consequence our 
simplified treatment is much more suitable for the latter because protons for 
sources within a radius ~ 50Mpc would be detected whereas photons would 
be mostly absorbed. 

It is clear that the feature of fig.2 is consequence of the hard component 
corresponding to quark hadronization reflected in the last term of eq.(14). 
If one should have considered an ordinary law 1/E^ at emission, the cor- 
responding flux due to uniformly distributed extragalactic sources would be 




(15) 
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given by eq.(14) without the last term and as shown in fig.3 would reasonably 
reproduce the observed events below the GZK cutoff but without recovery 
above it. 

We must note comparing figs. 1 and 2 that J is one order of magnitude 
larger for necklaces than for vortons which is consequence of the fact that 
the latter fit the flux at ~ lO^^e^ and the former that at ~ 10^^e\^ with 
a partial compensation due to degrading of energy. It is important to note 
that with the above values of r, /i and a monopole mass m ~ IQ^^GeV the 
energy per unit length due to monopoles turns out to be ~ lO'^'^GeV'^ and 
the fraction of critical density ~ 10~^, slightly smaller than that of ordinary 
strings. 



4 Conclusions 

We have seen that both vortons in halo and necklaces in universe may be a 
solution for the problem of UHECR above lO'^^eV. From the observational 
point of view the difference will be the expected anisotropy in the former 
case because of the asymmetric position of the sun in the galaxy compared 
to the isotropy of the latter characteristic of a cosmological origin. Regarding 
this, it must be noted that the anisotropy detected below the ankle is not 



observed above it|]TO| , which would be consistent with the appearance of 
an extragalactic component, a larger statistics being needed at the highest 
energy to see if new galactic sources contribute. 

Referring to elementary particle theory it is interesting that the most 
appropriate GUT models are different. For necklaces, since it is necesary 
that the breaking of the abelian symmetry leaves a discrete Z2 unbroken, a 
GUT model based on SO{10) is suitable. For vortons on the other hand the 
Eq gut model is better because the breaking of the contained additional 
abelian symmetry produces the necessary superconducting current with ex- 
otic fermions, whereas necklaces would not be formed due to the fact that 
its Higgs content does not allow an unbroken Z2 . 
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